Enhanced magnetic fluctuations in doped spin-Peierls systems: 

a single-chain model analysis 



Michele Fabrizio and Regis Melin 
International School for Advanced Studies (SISSA) Via Beirut 2~4, 34014 Trieste, Italy 

Abstract 

We analyze by means of real space Renormalization Group (RG) as well as by 
exact diagonalizations the properties of a single-chain model of a doped spin- 
Peierls system, where a major role is played by the localized moments created 
by the impurities. We are able to follow analytically the RG flow, which 
allows us to determine the relevant cross-over temperatures. In particular, we 
find an enhancement of magnetic correlations due to disorder, coexisting with 
an underlying dimerization, in an intermediate temperature range below the 
spin-Peierls critical temperature and above the coherence temperature of a 
regular array built by those localized moments (so-called soliton bandwidth). 
The possible relevance of these results to the doped inorganic spin-Peierls 
compound CuGeOs is discussed. 

I. INTRODUCTION 

The behavior of spin-Peierls systems in the presence of disorder has recently attracted 
considerable interest in the light of the intriguing properties of the inorganic CuGeC>3 com- 
pound. The pure compound is a quasi-one dimensional material, in which the interchain 
couplings in the two directions perpendicular to the chains are estimated to be 10% and 1% 
of the intrachain coupling. At at temperature T$p — 1AK, there is a structural transition, 
below which the CuC>2 chains dimerize and a gap in the spin excitation spectrum opens [Q. 



The properties below T$p are reasonably well described by a dimerized Heisenberg chain, 
with an additional next-nearest neighbour exchange coupling [0. 

If few percent of Cu is substituted by magnetic (Ni ||) or non magnetic (Zn [|],^]||) ions, 
or, if Ge is replaced by Si besides the structural transition, which still occurs close to 
14K, a second transition is detected at Tjv — 2 4- AK, which has been identified as a Neel 
antiferromagnetic transition. The staggered magnetic moment with 4% of Zn is estimated 
of the order of 0.2/ie [H- 

The simplest explanation of the appearence of a Neel phase upon such a weak doping, 
is that already the pure compound is quite close to a transition from a spin liquid phase, 
with a gap in the excitation spectrum, to an antiferromagnetic state, with gapless spin wave 
excitations. Then, one could imagine that doping effectively reduces the dimerization, or 
increases the interchain coupling, so that the system is pushed quickly in the magnetically 
ordered phase. Within this picture, one would represent the disordered system still as a 
homogeneous system but with modified parameters, and the enhancement of the magnetic 
fluctuations would be simply due to the critical quantum fluctuations. Notice that a dimer- 
ized exchange is not incompatible with a Neel order in more than one dimension, which 
explains the coexistence of dimerization and magnetic ordering. However, in our opinion, 
this explanation is not fully consistent with the experiments. A detailed discussion, including 
also a review of existing results, is postponed to the conclusion section. 

Therefore we believe that something else has to be invoked to explain why the spin- 
liquid phase of the pure CuGe0 3 is so unstable upon doping. In this paper we introduce 
and study another mechanism of enhancement of the antiferromagnetic fluctuations, which 
is purely due to disorder, exists also in a single chain, and can not be represented by a 
homogeneous system with modified parameters. Essentially, we will assume that, at low 
doping, the impurities release spin-1/2 solitonic excitations, which, however, are not free 
to move but get trapped in the vicinity of the impurities by interchain correlations P,|10 



Antiferromagnetic fluctuations can then be established since a coupling between these lo- 
calized spins is generated by the polarization of the spin-singlet background. However, one 



would naively disregard the effects of these magnetic fluctuations, since they are supposed 
to appear at an energy/temperature scale T* of the order of the spin- wave bandwith which 
would be obtained if these localized spins formed a regular lattice (soliton lattice band- 
width). This scale is proportional to the spin-Peierls gap with a prefactor exponentially 
small in the ratio of the average distance between the impurities to the soliton width. At 



2% Zn-doping, with the estimated soliton width of 13.6 lattice spacings |TT], this prefac- 
tor would be of order 0.025, while it would be 0.16 at 4% doping. Since at 2% Zn-doping 
T/v — 0.28Tsp ||, it would be hard to believe that these localized spins have something 
to do with the antiferromagnetism. We will show that this naive conclusion is not correct 
when these solitons are randomly distributed. In fact, as a consequence of disorder, sensible 
magnetic correlations will form well above Therefore, this effect may in principle play a 
role in the establishment of antiferromagnetism. Notice that this source of enhancement of 
magnetic fluctuations is not incompatible with the previously discussed scenario, rather it 
favours the approach to the critical point separating the spin liquid from the magnetically 
ordered phase. As we are going to discuss, we believe that there are evidences that these 
disorder-induced fluctuations play a role in these materials, especially at low doping. 

The paper is organized as follows. In section II we will introduce the model and discuss 
some existing results. In section III we will introduce the renormalization group approach 
which we use to study the model, and start the analysis, which will be further improved in 
section IV with the calculation of the spin-spin correlation function. In section V we present 
some numerical results which support the renormalization group analysis. Section IV is 
devoted to exact diagonalizations of small clusters, the interchain coupling being treated in 
mean field. The conclusions are given in section VII. 

II. THE MODEL 

The Hamiltonian which we consider in the absence of impurities is simply that of a 
dimerized Heisenberg chain 
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where 5 is the strength of the dimerization. As we said, to correctly describe the CuGe03 
compound, one should include a next nearest neigbour coupling. However, for our purposes, 
this is not really important. We assume that one impurity releases one spin- 1/2 solitonic 
excitation, connecting regions of different dimerization parity |12|]. The role of the interchain 
coupling is to provide a confining potential to the soliton, which will be trapped within some 
distance from the impurity |J (as well as to enlarge the soliton width [TB[]). Moreover, the 
weak link connecting the impurity nearest neighbors (which would be generated by virtual 
hoppings into and out from the impurity site, as well as by a next-nearest neighbor exchange) 
is approximated to be equal to the weak bonds in ([!]). This approximation is valid if one is 
interested really in what happens close to the middle of the spin-Peierls gap. Therefore, for 
a finite number n imp of randomly distributed impurities, the effective model will be assumed 
to consist of a squeezed chain with n imp sites less, described by the same Hamiltonian Eq.(fj]), 
but in the presence of randomly distributed domain walls. Their number will be < nj mp , since 
we can not exclude that pairs of solitons recombine. This amounts to take a site dependent 
5(i), which takes alternatively two values ±5, jumping from one to the other at the (random) 
position of the antiphase walls. If we assume no correlation between the impurities, then 
the distance r between two consecutive domain walls is distributed according to a Poisson 
law 

P(r) =ne- {r - a)n 6(r-a), (2) 

where n is the concentration of domain walls, and a is of the order of the lattice spacing. 
The XY version of this model [A = in Eq.([L])], has been recently studied in the continuum 
limit (a — > 0) in Ref. |TJ]]. In this reference, exact expressions of the density of states as 
well as of several thermodynamic quantities have been derived. These results are supposed 
to give a good qualitative description up to the isotropic point [0 < A < 1 in Eq. ([]])], since, 
as discussed in Refs. []T5J and [|nj, the spin anisotropy in a disordered chain does not play 
the same crucial role as in the absence of disorder. In this paper, we will study the low 



energy (lower that the spin-Peierls gap) behavior of the spin-isotropic model by means of 



the real space Renormalization Group (RG) originally introduced by Dasgupta and Ma |I7 



to cope with random spin chains, and further extended by Fisher ||16|| . We will show that 
this approach reproduces the exact results obtained in the XY limit in Ref. |[L4|| , thus proving 
the unimportance of the spin anisotropy for < A < 1. In addition, by this technique it is 
also possible to calculate the spin-spin correlation functions, which were not calculated in 
Ref. [|l4l ) being not easily accessible by the method used in that paper. 



A. Effective low energy model 

We start our analysis by building an effective low energy model which should describe 
the Hamiltonian (|l|) below the spin Peierls gap, and for low doping. In this regime, most 
of the spins are frozen into singlets apart from the spin-1/2 solitons localized around the 
domain walls. Two consecutive solitons are coupled by an exchange which we assume of the 
form 

J(r) = o e- r<fe , (3) 

where r is the distance between the two, which is distributed according to @, and l/<f>o is the 
soliton width. We have choosen units in which the energy has the dimension of an inverse 
length. Therefore the resulting low energy model describes a random Heisenberg model 
where spins are randomly distributed on a chain according to the law (|2|), and coupled 
by the exchange constants (0). This model should give a correct description of ([!]) at 
energies/temperatures smaller than O > an d for n <C 0o (the solitons should not overlap too 
much), and can be analysed by the RG technique of Ref. [|16j, as we are going to show in 



the following sections. A similar low energy picture has been proposed in Ref. |L8] on the 
basis of an exact diagonalization of a realistic model for CuGeOs, although limited to small 
clusters. By analyzing their data, the authors of that reference have concluded that the 
random Heisenberg model is indeed a good description for the low energy excitations. Here, 
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given for granted this low energy picture, we analyze some possible consequences, which 
can be useful for a comparison with experimental data as well as a guiding line for further 
numerical studies. 

It is now worthwhile to realize that the antiferromagnetic correlations in the squeezed 
chain do not simply translate into oscillations with wave vector Q = ir/a in the original 
chain. This is due to the presence of the impurities which change one sublattice into the 
other, since the coupling across the impurity site is antiferromagnetic. For instance, if in 
the squeezed chain a correlation function oscillates like cos(Qx), in the original chain it will 
oscillate like cos[Q(x — N(x))], where N(x) is the number of impurities in that interval. 
For randomly distributed impurities, after the averaging, this behavior will transform into 
cos(Qx)e~ 2ritX , where rij is the impurity concentration. Therefore, within this single chain 
model, we can not describe a true Neel long range order. For that, we would need to 
properly take into account the interchain coupling. A first attempt is discussed in section 
IV. C, where we calculate the susceptibility to a staggered (in the original chain, not in the 
squeezed one) magnetic field. However, for the purpose of describing how disorder enhances 
magnetic fluctuations, without pretending to push the analysis up to the AF transition, the 
single chain model we study is sufficient. 

To conclude, we notice that the XY version of this model is equivalent to a tight-binding 
Hamiltonian for spinless fermions at half filling with random hopping integrals. This model 
has been studied quite a lot in connection with the Anderson localization, for its intriguing 



properties. A detailed analysis was carried out by Eggarter and Riedinger [[L9| . They 
showed, quite generally, that the density of states has a Dyson singular behavior at the 
chemical potential pie) ~ l/|eln 3 e|, and that the localization length diverges logarithmically 
A(e) ~ | In e | . However, they did not discuss the consequences of these singularities in the 
correlation functions. This analysis was partially carried out by Gogolin and Mel'nikov 
20|j using the Berezinskii diagram technique. They showed, for instance, that, although 
one dimensional, this system has a finite conductivity. However, their approach is quite 
complicated and does not provide a simple interpretation of the low energy behavior. On 
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the contrary, the RG analysis carried out by Fisher for the spin model is, physically, more 
simple, and, once translated in the fermion language of the tight-binding model, provides 
not only a simple low energy picture but also many new results, as for instance the zero 
temperature power law decay of the average density- density correlation function. We have 
verified these predictions numerically in the tight-binding model, as we are going to discuss 
in section V. 

III. RENORMALIZATION OF THE BOND DISTRIBUTION 

We start giving a very short introduction to the RG transformation. Then, as a warm 
up exercise, we will calculate the renormalization of the probability distribution of the bond 
exchange- couplings . 

The RG transformation consists in successive eliminations of the most strongly coupled 
pairs of spins, by projecting out the Hilbert space onto the subspace where those pairs are 
frozen into singlets. The cut-off energy which is rescaled downwards in the course of the 
renormalization procedure is the strength of the strongest bond Q = max{J}. Therefore, 
going from Q to Q — SQ, amounts to project out all those pairs of spins which are coupled 
by an exchange of value Q. This operation generates an effective bond between two spins 
separated by a projected pair. If, for instance, at the energy scale Q, the coupling between 
spins 2 and 3 is J2,3 = then the projection onto the subspace in which those spins are 
frozen into a singlet generates a coupling between the neigbouring spins 1 and 4, given by 
Ji 5 4 = Ji,2^3,4/<^2,3- The decimation scheme becomes more and more valid in the course of 
the scaling procedure, as one verifies by the flowing of the bond probability distribution. 

In our particular previously discussed, we will consider a chain where spins 

are randomly distributed in such a way that the distribution of distances r between two 
consecutive spins is given by (|2[). The bond connecting two consecutive spins at distance r 
has a strength given by (|3]). Following Fisher, we define T = ln(0o/ max{J}), and a new 
variable ( = —T + ln((j) / J(r)). The recursion relation for the ^-variable is quite simple in 
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the above decimation scheme. In fact, in our previous example, where the sites 2 and 3 are 
decimated, a coupling between sites 1 and 4 is generated, such that £1,4 = £1,2 + C3,4 — C2,3 = 
Ci,2 + C3,4; being (2,3 = 0. We also define a scaling variable 77 = (/T. At our RG starting 
point, To = 0o<2, 77 = (r — a) /a, and the initial i] distribution is 

Q{v,^o) = ^pexp (-^77)0(77). 

With these notations, the RG scaling equation for Q is |L6 



r^ = Q + (i + 77)^ + Q(o,r) |VQ(?/,r)Q(77-7/,r). (4) 

We search for a solution of @ of the form: Q(?7, T) = /(T) exp (— f(T)r))9(r)), with / 
independent of 77. At the starting point /(r ) = na. The equation of motion of f(T) 



is 



-r^ = -/(r) + / a (r), (5) 

whose solution compatible with the boundary condition is 

0o (.1 ~~ TiaJ + 1 n 

If T = a0o, we recover the initial Q(j],To) distribution, and, in the fixed point limit, T — > 
+00, we have Q*{rj) = exp (—i])6(i]), which is exactly the random-singlet fixed point analyzed 
in Ref. In addition, we have an analytic expression of Q at all intermediate scales. 

If n(T) denotes the number of spins per unit length not yet decimated at scale T, we 
have p| 

dn(T) = 2Q(0,r) 

n(r) r 1 ' 

Using the previous solution for the bond distribution, taking the limit a — > 0, and denoting 
T = In (0 o /£) (E < O ), we have 

n(E) =n( - — ^- -) . (8) 

1 ; \nln^ /E)+<p ) 1 ; 



At the begining of the renormalization, E = </>q, and we have n(0o) — n, as expected. In 
the small E limit 

n(E) ~ 2 f° ; . (9) 

1 J nln 2 (</> /£) 1 ^ 

This leads to a low temperature uniform susceptibility 

y(T) = <H = (10) 

X{ ' AT 4nTln 2 (0 o /T)' 1 ' 



which is exactly the low temperature susceptibility obtained in Ref. [0 for the XY limit of a 
dimerized chain with a dilute random distribution of domain walls. The agreement between 
the two results is, first of all, a proof that the two models are equivalent at low energies 
(smaller than the spin-Peierls gap </>o), and, secondly, an evidence that the spin-anisotropy 
is not relevant. 

IV. BOND LENGTH DISTRIBUTION 

We now analyze the joint probability distribution of bond lengths and couplings, which 
allows to calculate the spin-spin correlations at arbitrary energy scale, and thus the relevant 
cross-over temperatures. 

Let P(C, I, r) be the joint probability distribution of a bond of length / and coupling ( 
at scale T. Following Fisher [|16|], the RG transformation of this quantity is 

dP((,l,T) d 



dr d( 

r+oo 



P(C,i,r) (ii) 



+ / dhdi 2 di 3 dCidCsS(C - Ci - CsW -h-h- h)P(o, h, r)P(Ci, h, r)P(C 3 , h, r). 

Jo 

For the Laplace transform of P((, I, T), 

P((,y,T)= dle~ l vp((,l,r), (12) 
Jo 

we look for a solution 

p(C,2/,r) = $( y ,r)e-^ r ^(c) (13) 



with $ and / two ^-independent functions. By ([□]), these two functions are solutions of the 
following differential equations 



dT 
dT 



$W), (14) 
/(y,r)<%,r). (is) 



Notice that (|13| ) is for the moment just an ansatz, which is valid only if one is able to solve 
([U]) and fll5|) consistently with the appropriate boundary conditions. For our particular 
model, the initial joint probability distribution is 

P(C, /, T ) = ne-^S (C - (I - a)0 o ) 0(C). (16) 

After Laplace transforming with respect to /, we find the following boundary conditions 
(a-0) 

/(y,r ) = / (y) = ^, (17) 

Tl 

$(y,r ) = -. (is) 

Solving Eqs.([L4]) and (|15|), we obtain 



r , [ru~\ M) coehtyC7(y)r) + ^(y) sinh(^(y)r) 
V / (y) Binh(Jd{y)T) + cosh(Tc^)r) 



$(y, r) = -JVc(y] _ (20) 

0o V / ( y ) sinh(^(yjr) + ^C^) cosh(^)r) 

where C(y) = f 2 (y, T ) — $ 2 (y, T ) = |/(y + 2n)/(j%. One can verify that the ansatz solution 
([jjD with the two functions defined by (p^)-([20|) does satisfy the appropriate boundary 
conditions. Moreover, when y = 0, we do recover the ^-probability distribution found in the 
previous section. If we scale y = y<pl/(2nT 2 ), and send r — > oo keeping y constant, we find 

K^H^ra e "*™> (21) 



which is the fixed point distribution found in Ref. [IB|. The joint probability distribution 
P(l, C, r) is obtained after the inverse Laplace transformation: 
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P(l, C, T) = / C+, °° ^<%, T)e-^9(C). (22) 

Having calculated P(l,£,T), we are now in position to derive various important quantities, 
as we are going to show. 

A. Fixed-point bond length probability distribution 

The bond length probability distribution P(l, T) at the fixed point is obtained by inte- 



grating (22) over £, and assuming y <C n. The result is |LC 



r) = ^ £ (n + i) (-1) V^V^, (23) 



where Z = l<^\j1n is the length in appropriate dimensionless units. This expression coincides 
with the probability distribution that a random walker remains inside an interval T after 
a "time" /. The same probability was found to play a crucial role in the study of the 



one-dimensional tight-binding Hamiltonian with off-diagonal disorder |19 |. The comparison 
between the RG approach and the analysis of the Schrcedinger equation for the tight-binding 
Hamiltonian provides the quantum-mechanical interpretation of the bond length probability 
distribution at scale T as the probability distribution of half of the distance between two 
consecutive nodes of the wavefunction at energy E = 0oexp (— T). The average bond length 
or, equivalently, half the average distance between the nodes at energy E, is simply given 
by 

2nln 2 (0 o /£) 



< I >- 



from which one obtains an estimate of the integrated density of state N(E) = 1/(2 < / >), in 
perfect agreement with the exact result obtained for the dimerized XY chain with randomly 



distributed domain walls 14 . 



B. Spin-spin correlation functions 

The spin-spin correlation function 
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Xa (l) = (S a (0)S a (l)), 

(a = x, y, z), can be also calculated by making use of the joint probability distribution. At 
finite temperature T < <fi , the RG has to be stopped at an energy scale IV = ln(0 o /^). 
If a bond of lentgh I has a coupling larger than T, it will be decimated, and contribute a 
constant to x(0- n the contrary, for the bonds of the same length but coupling smaller 
than T, one can, as a first approximation, perform a high temperature expansion, keeping 
only the first non vanishing term. As a result, we assume the following expression for the 
spin-spin correlation function: 

X (i,T)~ J d(^p-n(T T )P((,l,T T ) + dTn(T)P(0,l,T) 

= n(r T ) J d(e-<P((,l,r T )+ dTn(T)P(0,l,T) = xi(l,T) + X 2(l,T), (24) 

where xi is the result of the high temperature expansion (notice that, by definition, £ = 
ln(T/ J) at scale IV), and xi is the contribution of the decimated bonds. Notice that, by the 
way in which the decimation scheme is built, the distance of a bond connecting two spins 
is always odd in units of the lattice spacing, at least in the squeezed chain. Therefore the 
above correlation function is in fact staggered in the squeezed chain. 

Since a very long bond will most likely form at scale r < IV 10, Xi wm dominate 
the correlation function when I is very large and for intermediate IV. (On the contrary 
for IV — ► oo it is X2 which dominates.) In this limit, x — Xi slightly differs from the 
expression used by Fisher, i.e. x(0 ~ P{^J^t)- I n fact, the difference is irrelevant at low 
temperature, but is important at higher temperature. In particular, with the definition (p4|). 
we can recover at T ~ </> the correct asymptotic behavior x(Z) ~ exp (—<f) l), for O 3> n. 
Since we are mostly interested in the cross-over, we will use fl24]). 

1. Thermal correlation length 
Let us first discuss the behavior of Xi i n (|2~Sp. After the ^-integration, one has to calculate 
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Xi{l,T)=n{T T ) v 25 

Jc-ioo 2m 1 + j{y, Ft) 

At finite temperature, the long distance behavior will be dominated by the nearest singu- 
larity to the origin in the complex y-plane. This pole will define the inverse of the thermal 
correlation length £r, by ~ e~ l ^ T . We will compare the resulting £t with the thermal 
correlation length £ t of a chain where spins are regularly distributed with interspin distance 
\jn and coupled by an exchange Jo = (poe~^°^ n (soliton lattice). For convenience, we will 
consider the longitudinal spin-spin correlation function in the XY limit of both models. This 
amounts to substitute 1 with 2 in the denominator of Eq . (p5|) . This comparison is useful to 
understand in which regime the spin-spin correlation is reduced by disorder, and in which 
regime, if any, is on the contrary enhanced. We find that 

?„r = -, - . (26) 

2«l"(fW(S) + 1 ) 

If T <C T*, where T* = 2 Jo is the bandwidth of the spin excitations, then 



2rdn 2 (0 o /T) 

^ = ~~^o ' ( } 



The correlation length diverges as T — > 0, even though more slowly than for the regular 
lattice, where £ot ~ Jo/(imT). Therefore, for very low temperature, we find that the 
disordered system has weaker correlations than the ordered soliton lattice, a result rather 
obvious. 

On the contrary, within the interval T* < T < O ? a different situation is encountered, in 
which the correlation length of the disordered system is larger that the one of the ordered 
lattice, as shown in Fig. |l[ We find that the correlation length of the disordered system 
increases linearly in ln(0 o /T) in that temperature range, while £o remains almost constant. 
This result is also somewhat predictable, even though at first glance surprising. In fact, 
above the coherence temperature T*, the disordered system is more correlated since it takes 
advantage from configurations where the spins are closer than in an ordered lattice. In 
connection with our original model of a doped spin-Peierls system, this result implies that 
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magnetic correlations may appear at temperatures well above the soliton bandwidth, which 
is not in contradiction with the experimental evidences in CuGeOs. 

Let us now consider the limit T — > , where the contribution of \2 in ( F^l) is dominant. 
With exponential accuracy, we can use the fixed-point joint probability distribution (pT|). 
After inverse Laplace transform, 



2 J. oo 



p(o,/,r) = ^ E(-i) m+1 ^ 2 ^ 2m2[/r2 - 

Performing the integral over T and then summing the series, we find x(0 — n /{3<pol 2 ), which 



is the known power law behavior [16] for a random Heisenberg chain. 



C. Staggered magnetic susceptibility 

Another relevant quantity which can be calculated is the staggered magnetic suscepti- 
bility. We add to the Hamiltonian a term IiJ2r(— 1) R S z r, where R is the localized spin 
position in the original chain, and calculate the second order correction to the free energy 

5F(T) = -— / dn / dr 2 £ (~l) R+R « S zR {t x )S zR ,{t 2 ) », (28) 
L Jo Jo RRI 

where << • • • >> denotes a thermal and disorder average, and L is the size of the chain. 
At scale T, mostly the spins which are still free will contribute to (p8|). Therefore, one can 
again perform a high temperature expansion for these spins and keep only the zeroth and 
first order term. 

At zeroth order, only R = R' contributes to (p8[), and gives 

which results into a staggered susceptibility equal to the uniform susceptibility of Eq. (|i~0|). 
This clearly reflects the fact that only on-site correlations are involved, which do not distin- 
guish between uniform or staggered fields. 

Longer range correlations start to appear at first order. We find 
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>, 



where the remaining average is over the disorder, and Jr,r>(T) is the exchange between the 
spins at sites R and R' in the original chain at scale T, on provision that these two spins 
have not been decimated yet. At this point we have to be very careful in distinguishing the 
original from the squeezed chain. Notice that the number of spins which have been already 
decimated and which lye in between the two spins is by construction even. If we assume 
that each spin is localized very close to each impurity, then also the number of impurities 
between R and R' will be even. Therefore the parity of the distance R — R' between the two 
spins depends whether they are at the right or the left of the impurities to which they are 
bound. In particular, the case in which the distance between the two spins in the squeezed 
chain is I can correspond to four possible cases in the original chain. The first two cases, 
which correspond to R — R' even, are when R and R' are both on the right or on the left 
of the impurities to which they are bound. We define these two cases, in obvious notations, 
as (+, +) and (— , — ). The other two cases, corresponding to R — R' odd, occur when one 
spin is on the left and the other on the right and viceversa, which we define as (— , +) and 
(+, — ). The value of the coupling will be different in the four cases. What enters in the 
calculation of the staggered susceptibility is the combination J_ ;+ + J +; _ — J_ _ — J +i+ of 
the exchange constants in the various cases. Due to the exponential dependence of J upon 
the distance, this combination will be always positive, as if R — R' were effectively odd. We 
therefore assume that 



\R+R' 



Jr,r>(T) = 1 J 1 (T)>0, 



where Ji(T) is the exchange coupling of a bond of length I in the squeezed chain at scale T, 
and < 7 < 1 is a reduction factor. With this assumption we get 

iFm[ T )=7 ^T ) J^P ( UT T )=^ f J^ } , (29) 
where / has been defined in Eq.(|6]). At low temperature, this gives a contribution to the 
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staggered susceptibility ~ 1 / (T"| lnT| 3 ), which is subleading with respect to the zeroth order 
term, but still diverging. 

V. NUMERICAL CALCULATIONS AT ZERO TEMPERATURE 

We now turn to the numerical calculations of the zero temperature correlation length of 
our effective model. In the framework of the renormalization group approach to the bond 



disordered model, Fisher [16[ showed that the zero temperature correlations << S"Sj >> 
behave like — j| 2 , whatever a = x,y,z, as we have seen in the previous section. Our 
goal is to check this prediction by means of exact diagonalizations in the XY case. 



We use the Jordan- Wigner transformation |2T] to map the spin 1/2 problem to a spinless 



fermion problem. We work with an even number of sites so that the fermions are periodic, 
and in the S z = sector, implying half-filling. The XY Hamiltonian in the fermion language 
is the one of non interacting fermions with random hoppings 

H = \H^{4 + iCi + cfc l+1 ), (30) 

i 

where the Jj exchanges are all antiferromagnetic and distributed at random. We choose 
the bond distribution P(J) = @(J)@(1 — J). For this model, we will calculate both the 
<< S*Sj » and << S^~Sj~ » correlations. This model belongs to the same universality 
class as the model defined by @ and @. In order to test the equivalence between these 
two models, we will calculate the ((SqS^)) correlations of the model (fj) and (§) and show 
that it indeed decays like 1/R 2 . 

A. S Z S Z correlations 

The first step is to calculate numerically the spectrum and the eigenvalues of the tight- 
binding Hamiltonian 

# = ~£Ji(|i + iM + |i)(z + i|). (3i) 
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Let |\l/ a ) = Ya the eigenstates of (|31~]). The quantum average over the half-filled Fermi 
sea (FS) is calculated as < S Z S Z >= A iyj - Bi/2 - Bj/2, with 

■ 1,, E (^?) 2 +E E (32) 



a,/3eF5 a£FS f3£FS 

Bi = E (*?)' • (33) 

The results for the << SfiSJ >> correlation are shown on Fig. 0, in good agreement 
with the — j\ 2 scaling. 

A main feature of the random singlet fixed point is that typical and average correlations 
differ. Fisher [16] claims that the typical correlation behaves like — ln|C*J p | ~ \i — j\ 1 ^ 2 . 
In order to check this prediction numerically, we calculated the typical S Z S Z correlation 
function, (defined as the most probable correlation) as a function of distance. The quantity 
— In \Clj P \/\i — jl 1 / 2 is plotted on figure [3], where a plateau is clearly visible, indicating a 
consistency between Fisher's prediction and our numerical calculations. 

In the case of the (0) and (H) model, the << S Z X S Z X+R » correlations are shown on 
figure [|, where a 1/R 2 behavior is also visible at large distances. 



B. S + S correlations 



We now turn to the computation of the S + S~ correlations of the XY chain. Following 
we define Ai = cf + q and Bi = cf — c\. Assuming % < j, expanding the string 
operator in terms of the A and B operators and performing the adequate contractions, we 
get < SfSj >= Dij/2, with the determinant 



Da 



Gi,i+l Gi,i+2 ••• Gi,j 

Gj-i,i+\ Gj-i t i + 2 ■■■ Gj-ij 



(34) 



where 
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G id =< BiAj >= E " E (35) 

\a&FS a^FS/ 

This allows us to calculate numerically the << S^Sj~ » correlations of the disordered XY 
chain. The results are plotted on figure [5|, where a cross-over to a 1/R 2 behavior is visible. 
Notice that in the S + S~ case, the 1/R 2 behavior is achieved for larger separations than in 
the S Z S Z case. In addition, we have also calculated the equal-time single particle Green 
function, which also seems to decay as 1/R 2 at large distances. 



VI. DIAGONALIZATION OF SMALL CLUSTERS 

In order to have a further confirmation of the physical picture that comes out from the 
previous analysis, we have performed exact diagonalizations of small clusters. 

We work on a chain of 12 sites with 2 non magnetic impurities. We enumerate the possible 
disorder realizations, which are shown on figure || and diagonalize the spin Hamiltonian in 
all the S z sectors. We also make the assumption that the domain wall created by each 
impurity in the dimerization pattern does not move, but remains trapped at the impurity 
site. This simply implies that each impurity releases a free spin which is located either at 
the right or at the left of the impurity site (see the figure |^). A strong bond has a value 
J(l + 5) and a weak bond J(l — 5). Since we work with a very small chain, the dimerization 
should be large enough to make the spin-Peierls correlation length much shorter than the 
average spacing between impurities. The Hamiltonian parameters are thus far from being 
realistic. Moreover, we assume that the second nearest neighbor interaction across the 
impurity is equal to the weakest bond J(l — 5). In practice, we take J = 1 and 5 = .6. 
The reason why we cannot address the problem for realistic values of the parameters is that 
realistic computations would involve a much larger number of spins, which is not accessible 
numerically. However, as we are going to show, even our over-simplified system leads to 
some interesting results, not in disagreement with available experimental data. 

In order to describe the antiferromagnetic transition we treat the interchain coupling J± 
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in a mean field approximation. Specifically, we impose a staggered magnetic field h s and 
calculate the staggered magnetization in the presence of h s . We thus get the values of the 
staggered magnetization ms(hs,T) as a function of the staggered field and temperature T. 
Next, we determine h s {T) self-consistently by imposing that hs(T) = 2J±ms(hs,T). The 



From that figure we see that the on-set temperature at which a finite staggered magnetiza- 
tion appears (in a sense the mean field Neel temperature) is different for different disorder 
realisations. In particular, while most of the realisations have almost the same on-set tem- 
perature, there is one which has a larger temperature. This confirms our previous result 
that magnetic correlations start to appear above the soliton bandwidth and are due to rare 
disorder realisations which however dominate in that temperature range. 

Next, for each disorder realization, we have calculated the susceptibility at the self- 
consistent point and averaged over all the disordered realizations. The average susceptibility 
as a function of temperature is plotted on figure |[ As the temperature is decreased, the 
susceptibility starts to decrease for temperatures smaller than the spin-Peierls gap. Instead 
of going monotonically to zero for decreasing temperatures, as in the case of a pure spin- 
Peierls system, the susceptibility reaches a minimum, after which it start increasing, and 
finally is cut-off by the antiferromagnetic order. This behavior is in qualitative agreement 
with the experimental data [0. 

We also plotted on figure |^ the variations of the average overlap with the antiferromag- 
netic state. If |\?af) denotes this state, then the overlap is defined by 



where a labels the eigenstates |\& a ), and E a are the energies. As expected, we observe a 
sharp increase of this overlap as the self consistent staggered magnetic field is switched on. 
However, we also observe an increase, although slow, at larger temperatures. This increase 
is consistent with the existence of enhanced antiferromagnetic fluctuations above the Neel 
temperature that were described in this paper. 



self-consistent staggered field is plotted on figure |7] for the different disorder realisations. 




(36) 
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VII. DISCUSSION AND CONCLUSIONS 



In this article we have analyzed both analitycally and numerically a model for a doped 
spin-Peierls system. In this model a major role is played by the localized spins which are 
released by the impurities. These spins are antiferromagnetically coupled by the polarization 
of the singlet background. An important result that we find is that sensible antiferromagnetic 
fluctuations start to appear just below the spin-Peierls gap and well above the average value 
of the exchange coupling, which is related to the coherence temperature of an ordered lattice 
made by those spins. In other words, there is a high temperature regime where the disordered 
system is more correlated than an ordered one. In this temperature range, we find that the 
thermal correlation length increases logarithmically with temperature. Moreover, in spite 
of the presence of disorder, these fluctuations are quite long ranged and, in fact, give rise 
at zero temperature to power low decaying spin-spin correlation functions, which has been 
numerically verified in the XY limit. 

We believe that these disorder-induced magnetic fluctuations are important in doped 
CuGeOs, especially in driving the system towards a Neel ordered phase. There are exper- 
imental evidences which, in our opinion, reveal the importance of these fluctuations. For 
instance, the data of Ref. || show no saturation of the intensity of the antiferromagnetic 
Bragg peak as the temperature is lowered below T/v- This has been interpreted by the au- 
thors as a clear signature of disorder. In addition, the extrapolation at T = of the value 
of staggered magnetization in the Neel state, which is 0.2//g for 4% Zn-doping according to 
Refs. II and ||, implies that not all the Cu spins are involved in the antiferromagnetism. 
Finally, we find interesting the uniform magnetization data M(T) of Ref. |7j], for Si doped 
compounds. In fact, for doping < 0.005, with decreasing temperature, M(T) shows a first 
drop at Tsp, followed by a minimum, after which M(T) starts increasing, without showing 
any saturation at low temperature. This is a clear evidence that doping induces some lo- 
calized moments which are randomly distributed. For larger doping, the low temperature 
singularity is cut-off by a second transition into a Neel phase. These results suggest that 
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those localized moments do participate actively to the antiferromagnetic ordering. It is 
moreover suggestive in that Reference that the doping concentration at which Tjy is maxi- 
mum seems to coincide with the vanishing of Tsp, or better the closing of the spin-Peierls 
gap, as if, in the presence of an underlying spin-gap, doping helps antiferromagnetism, which 
is instead reduced by doping in the absence of the gap. This result is compatible with our 
model. 

Unfortunately, we can not give a realistic estimate of the Neel temperature within our 
single-chain toy model. This would need the inclusion of interchain couplings which can not 
be handled by the RG approach we have used. 

It is now worthwhile to discuss more in detail our theoretical motivations in support of 
these disorder-induced magnetic fluctuations. As we said in the introduction, an alternative 
explanation of the appearence, upon such a weak doping, of a magnetically ordered phase, 
is that already the pure compound is quite close to the quantum critical point separating 
the spin-liquid phase, with a spin-gap in the excitation spectrum, from the Neel phase, with 
gapless spin wave excitations. Then one can imagine that doping effectively decreases the 
dimerization or increases the interchain coupling so to push the system in the Neel phase. 
However, this simple view is in our opinion not fully satisfactory, as we are going to discuss. 

Katoh and Imada p3| have recently analyzed by quantum Monte Carlo and spin wave 
theory a quasi one-dimensional Heisenberg model with a dimerized exchange. Their model 
differs from the one expected to describe the CuGe03 for the absence of an intra-chain next- 
nearest neighbour coupling, as well as for the fact that in this compound the dimerization 
is staggered in the transverse direction, contrary to the model analyzed by those authors. 
Anyway, what Katoh and Imada find is that, for a realistic value of the ratio of the interchain 
to intrachain exchange (~ 0.1), the critical value of the dimerization separating the Neel 
from the spin-liquid phase is 5 C ~ 0.05, not far from the estimated value of 5 = 0.03 in 
the pure CuGeO"3. What is more important for what we are going to say, is that spin- wave 
theory overestimates 5 C , giving a value ~ 0.7. If one does spin-wave theory for a realistic 
model, taking into account that the dimerization is staggered in the transverse direction 
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and including a next-nearest neighbour coupling, one finds that a model describing CuGeC>3 
is well inside the spin-liquid phase, and it would need a really large modification of the 
Hamiltonian parameters to enter the Neel phase. One could imagine that the quantum 
fluctuations not present in the simple spin-wave theory might change these results and push 
the system towards the Neel phase. However these quantum fluctuation corrections go in 
the opposite direction in the case analyzed by Katoh and Imada, and we do not understand 
why they should work oppositely in the realistic case. 

There have been other attempts, which make use of bosonization, to derive the phase 
diagram of a quasi-one dimensional Heisenberg model with dimerization and intra-chain 



next-nearest neighbour exchange |24],13[|. In Ref. it is claimed that, within this approach 
it is possible to show that the model describing the pure CuGe0 3 is indeed very close to the 
transition towards a Neel phase. However, also these approaches are not fully satisfactory in 
our opinion. First of all, as a consequence of the one- dimensional treatment (the interchain 
coupling is treated in mean-field, assuming that the neighbouring chains provide a staggered 
magnetic field, which is then calculated self-consistently), the transition between the spin- 
liquid phase and the Neel phase is accompanied by the vanishing of the dimerization. This is 
not what is seen experimentally (see e.g. Ref. 0). In fact, as we already said, a Neel phase 
is not incompatible with a staggered exchange in more than one dimension. A simple spin- 
wave calculation shows that there is a gain in energy by switching on a small dimerization 
S, given by 

1 



E(6)-E(0) = ^Ev(' / + J !/- 2 J' sin 2 k x ) 2 - ( J cos k x + J y cos k y ) 2 - J 2 5 2 sin 2 k x 
k 



- J (J + J y - 2 J' sin k x ) 2 — (J cos k x + J y cos k 



) 2 

y) ' 



where J(1±S) is the intra-chain nearest neighbour exchange for bonds of different parity, J y 
is the interchain exchange (we have assumed a two-dimensional system, since the coupling 
in the other direction is much smaller than J and J y ), and J' in the intra-chain next nearest 
neighbour exchange. This gain is quadratic in 5, for small 5, and should be compared with 
the loss of elastic energy due to the lattice distortion, which is also quadratic. Therefore, it 
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is not impossible that a Neel state is stable in spite of a dimerized exchange. In addition, 
the bosonization approach is extremely delicate and not fully justifiable for a J' ~ 0.23J, 
very close to the value 0.2411 which separates the gapless Luttinger liquid phase from the 
gapped valence bond regime. 

In conclusion, we believe that one really needs to take into account those disorder induced 
magnetic fluctuations to explain why the spin-liquid phase of the spin-Peierls compound 
CuGe0 3 is so unstable upon doping, and easily gives way to a Neel phase. The analysis 
presented here is a first attempt in this direction. 
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FIG. 1. Thermal correlation lengths of the disordered system with n/tfto = 0.1 (solid line), and 
for the non disordered one (dashed line). 



26 




FIG. 2. Average correlation for N = 26, TV = 50 and N = 102 sites. 200000, 350000 and 
120000 disorder realizations were used respectively. The solid line is the fit to the 1/x 2 behavior. 
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FIG. 3. Behavior of the typical correlation as a function of the separation. We have plotted 
— In \C*j P \/\i— j'l 1 / 2 as a function of the spin separation, calculated for a N = 102 chain and 120000 
disorder realizations. A plateau is visible for distances between 20 and 35, consistent with Fisher's 
prediction for the large scale behavior of the typical correlations. 
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FIG. 4. S Z S Z correlations of the model @ and © for = 1 and 1/n = 5, and iV = 12.JV = 26, 
iV = 50 and iV = 102 sites. 50000 disorder realizations were used excepted for N = 102 where 
13000 disorder realizations were used. 
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FIG. 5. S+S~ correlations of the disordered chain for N = 50, iV = 102, N = 202, N = 402 
chain. The number of disorder realizations is 1000 except for the N = 402 case where 420 disorder 
realizations have been used. The solid line indicates the l/R 2 behavior. 
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FIG. 6. The eight inequivalent disorder realizations of the 12 sites chain with 2 impurities. The 
empty circles are the impurities, the filled symbols are the spin 1/2 sites, the filled squares being 
the solitons. A double line stands for a strong bond (J + 5) and a single line stands for a weak 
bond (J - 5). 
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FIG. 7. Staggered magnetic field as a funtion of temperature for the eight disorder realizations 
(J_l = 0.315). 
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FIG. 8. Susceptibility versus temperature (J± = 0.315). 
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FIG. 9. Variations of the overlap with the antiferromagnetic state as a function of temperature 
. = 0.315). 
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